A quantitative evaluation of metallic conduction in conjugated polymers 
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As the periodicity in crystalline materials creates the optimal condition for electronic derealiza- 
tion, one might expect that in partially crystalline conjugated polymers derealization is impeded 
by intergrain transport. However, for the best conducting polymers this presumption fails. Dereal- 
ization is obstructed by interchain rather than intergrain charge transfer and we propose a model of 
weakly coupled disordered chains to describe the physics near the metal-insulator transition. Our 
quantitative calculations match the outcome of recent broad-band optical experiments and provide 
a consistent explanation of metallic conduction in polymers. 

PACS numbers: 71.23.An, 71.30.+h, 72.15.Rn, 72.80.Le 



One of the hallmarks of metallic transport is the nega- 
tive dielectric function e at low frequency u>, that reflects 
a time lag between induced current and applied field due 
to the inertia of delocalized charge carriers. The surpris- 
ing observation of negative e in conjugated polymers by 
Kohlman and co-workers^ which was recently verified by 
Romijn et al. and Lee and Heeger spurred experimen- 
tal and theoretical research into the nature of metallic 
conduction in these materials. As displayed in Fig. ^ 
the full spectral response is quite complicated including 
multiple zero crossings of s(u>) and non-monotonous be- 
havior of the conductivity cr(uj). The low-w dynamics are 
characterized by long scattering times r ~ ps and low 
plasma frequencies lu p ~ meV. This differs three orders 
of magnitude from conventional metals where r ~ fs and 
w p ~ eV, while the carrier density n is only one order of 
magnitude less. An empirical scaling relation between r 
and w p has been found, see Fig. |3 suggesting a common 
mechanism governs these parameters. 4 
The initial models for metallic polymers relied either on 
the Anderson theory of localization^^ or on a granular 



picture,^ but both predict e 
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Including per- 



colation effects explains negative but the calculated 
Up lies two orders of magnitude above the experimental 
values, 9 Prigodin and Epstein suggested that the metallic 
state is sustained by atypical resonant tunnelling events 
and used the Landauer-Buttiker transmission framework 
to explain small lo p and long r. 10 However, the conditions 
for resonant tunnelling are not always fulfilled^ and the 
detailed physical mechanisms remain to be untangled. 
Recently, Prigodin and co-workers successfully explained 
the w-dependence in the insulating phase using quasi- ID 
variable-range hopping theory with interchain transfer as 
rate-limiting step^i By accounting for quasi-lD conduc- 
tion, Kaiser et alH gave a description for metallic con- 
ductivity of granular polymers, which works for blends as 
welli 1 ^ Recent optical experiments by Lee and Heeger^ 
and us^ could be explained in a ID picture as well. Also 
Kohlman et al. discussed ID effects^ Here, we attempt 
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FIG. 1: Optical dielectric function and conductivity of metal- 
lic poly pyrrole. Symbols: experimental results taken from 
Ref. 2. Drawn lines: theoretical prediction for weakly cou- 
pled disordered chains (see text). 



to elucidate the nature of metallic conduction in polymers 
in more detail. We first argue why the metallic phase is 
dominated by the ID structure. Then, like Ref. ITot we 
use a Landauer-Biittiker approach to model the transport 
and show how interchain coupling drives the intrinsically 
insulating ID phase to a metallic state at the cross-over 
from ID to 3D coherent conduction. All transport pa- 
rameters can be described in terms of disorder and cou- 
pling strength and the relevant length scales are £ the 
ID localization length, and A the length over which car- 
riers experience sequential interchain events. For weak 
coupling £ < A, transport keeps a ID signature and ef- 
fectively only a fraction of the carriers is involved in 3D 
transport. Our model reproduces the reported experi- 
mental results well and naturally ascribes the unusual 
low-w dynamics to small interchain coupling. 

In metallic conductors, delocalized states exist at the 
Fermi level. Delocalization is counteracted by disorder 
and the balance between metallic and insulating ground- 
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FIG. 2: Scaling relation between plasma frequency and 3D 
scattering time. Dots: experimental results taken from Ref. 4. 
Full line: calculation for disorder driven MIT (see text). 



states is critically dependent on the weakest charge trans- 
fer steps in the system. In the disordered polymers, 
chains are aligned within microcrystallites that are sep- 
arated by amorphous regions, and derealization may be 
impeded by either intergrain charge transfer or by inter- 
chain charge transfer. To identify the critical mechanism, 
we compare the energy scales. 

Intergrain charge transfer - For the best conducting 
polymers the crystalline fraction is f c — 50% (polypyr- 
role, PPy, and polyanilinc, PAni), even up to f c = 90% 
for polyacetylene (PAc); the typcial microcrystallite size 
D c amounts to 2 nm (PPy), 5 nm (PAni), or 10 nm 
(PAc). For isolated grains^ geometrical considerations 
imply a width of amorphous regions W a — D c (f c — 1), 
giving W a — 0.5 nm, 1.3 nm, and 0.4 nm, respectively. 
Prigodin and Epstein estimate the localization length 
along the chain ( = 1.2 nm>i£ A different approach is 
based on the potential barrier height B of the amorphous 
region, which from transport activation energy measure- 
ments we estimate B = 0.01-0.1 eVii* A wavefunc- 
tion penetrating this barrier falls off exponentially with 
C = (2m c B/h 2 )- 1 / 2 = 0.6-1.9 nm. This corresponds to 
fully back-scattered waves, expected for strong disorder, 
and agrees with the experimental value £ = 0.8 nm for 
insulating polythiophcne, 15 Another approach is based 
on a result by Thouless: £ = a(Ep/B) 2 ^ which easily 
yields 10 nm using reasonable parameters (Ep = 0.5 eV, 
B = 0.1 eV, and a — 4 A the microscopic length scale, 
e.g. a monomer). Taking a conservative estimate £ = 
1.2 nm, the intergrain transfer integral through amor- 
phous regions along ID chains I\\ = Iq exp(— 2W a /Q), 
with Iq « 2.5 eV the intra-chain 7r-electron transfer of 
the unperturbed chain, equals at least 1.1 eV, 0.3 eV, 
and 1.3 eV for PPy, PAni, and PAc respectively. 

Interchain charge transfer - In conjugated polymers, 
the 7r-electron system is formed by hybridization of 2p x 
orbitals on adjacent carbon atoms that constitute the 
polymer backbone. Interchain charge transfer stems 
from the overlap of the 7r-electron clouds on neighbor- 
ing chains. We estimated the interchain charge trans- 



fer I± for parallel chains of carbon 2p x orbitals using 
the Slater wavefunction approximation. Our calcula- 
tions suggest an optimum packing distance of 3 A where 
7j_ ~ 10 meV is maximized. For larger separation, I± 
falls off exponentially with interchain localization length 
~ l A in agreement with experimental results^ Ex- 
tensive tight-binding calculations by Mizes and Conwell 
yielded I x ~ 30 meV for PAcii 

Thus, for f c « 50 % the energy scale for interchain 
charge transfer is at least an order of magnitude smaller 
than that for intergrain charge transfer^ This is not 
surprising when it is realized that intergrain transport 
is essentially an intrachain process that proceeds along 
the covalently bound chains, while adjacent chains are 
not chemically but Van der Waals-bonded and hence in- 
terchain electronic coupling should be much weaker. It 
means that (i) intrachain carrier derealization extends 
over several grains and (ii) formation of truly delocalized 
states is governed by the small interchain charge trans- 
fer. Therefor, we propose a system of weakly coupled 
disordered chains to model conducting polymers. Our 
approach is as follows: we limit ourselves to an isolated 
disordered ID chain where quantum mechanics leads to 
localized ID wavefunctions ^id- We argue that, stay- 
ing within a ID framework, interchain carrier exchange 
is equivalent to dephasing of the ^ids and derive the 
condition for derealization. Finally, we apply our model 
to describe the carrier dynamics of conducting polymers 
using known microscopic parameters and a tight-binding 
approximation for the interchain coupling. 

First, we briefly review the case for isolated chains. 
Even for minimal disorder, the electronic eigenstates ^id 
are localized due to repeated coherent back-scattering.- 
A convenient approach to solve the transport problem 
is the Landauer-Biittiker transmission framework^ As- 
sume a transmission coefficient T < 1 per unit length a. 
Then, over a distance z > a the envelope of an unreflected 
wave falls off as T z l a implying an intrachain localization 
length Q = a/ln(l/T). Disorder determines the micro- 
scopic transmission probabilities and hence sets C- The 
mean level spacing of the localized states is A = l/(g() 
with g the density of states per unit length and unit en- 
ergy. The dimensionless conductance G = T/(l — T) and 
it follows that at short length scales (z < £) G is finite 
and given by the classical (Ohmic) result G(z) — (/z, 
while for a long chain G(z) = exp(— z/£) scaling is 
non-Ohmic and the conductivity, Gz, becomes zero for 
z — > oo. The conductivity due to resonant transitions 
between localized states is appreciable only for hw ~ A. 
For hw > A, this Lorentz oscillator response approaches 
the Drude conductivity with plasma frequency Q p ~ eV 
and intrachain scattering time tq ~ fs. 

The above discussion is essentially based on a one- 
electron Schrodinger equation assuming phase-coherent 
transport along the chain. Interchain transfer introduces 
a probability e for interchain transitions and formally the 
isolated chain eigenfunctions are no longer appropriate. 
For small e, however, the isolated chain picture remains 
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accurate for short time and length scales: the \tiDS ob- 
tain a finite lifetime. In case of strong coupling such a 
picture is not realistic and the carrier transport prob- 
lem should be considered for the fully coherent 3D sys- 
tem. In conjugated polymers interchain charge transfer 
is weak and a ID transport model including the effect 
of coupling seems a good starting point. Particle con- 
servation implies that interchain coupling occurs in the 
form of carrier exchange with a reservoir, where the car- 
rier's initial state \Pid.i differs from the final state 'J'id^, 
which thus introduces the finite lifetime of the ID states: 
interchain transport leads to dephasing reflecting that 
the 4"ids are not "true eigenstates" of a system of cou- 
pled chains. Usually, phase-breaking is associated with 
inelastic events occurring at non-zero temperature. How- 
ever, phase-breaking is not necessarily dissipative when 
the initial and final states have the same energy^ Al- 
lowing "elastic phase-breaking" implies that the "true 
eigenstates" are extended over the entire system of cou- 
pled chains since then the energy level spacing b 2 /(gL 3 ) 
vanishes for L — > oo (6 interchain distance), i.e. the 
system is gapless, or in other words a metal. The condi- 
tion for which this becomes possible determines the zero- 
temperature MIT and will be discussed further below. 

Phase-breaking events are modelled in the Landauer- 
Biittiker framework by incorporating current-conserving 
phase-randomizing scatterers into the chain: each scat- 
terer randomizes phase with amplitude y/e yielding prob- 
ability e for phase-breaking. The total transmission prob- 
ability is given by T tot = T c + T; ^ with T c the direct co- 
herent intrachain transmission probability (no interchain 
events) and T; transmission including at least one inter- 
chain transition. The conductance 

G = Ttot = - + - (1) 

1 - T tot 1 - T c - T; 1 - T c - Ti V ; 

consists of two contributions: proportional to T c and T;, 
respectively. However, the relations are not straightfor- 
ward, as the denominators contain both T c and Ij. 

The carrier fraction e that suffered dephasing at z = 
will, on average, suffer a subsequent dephasing event at 
z ~ A = a/e. We can discern two limits: z < A and z > A. 
When z < A, T c > T- - 0(e), Eq. © simplifies to that 
of the isolated chain: G = T c /(1 — T c ). When z > A, T c is 
exponentially small and G — Ti/ (1— Ti). Strong coupling, 
i.e. every scatterer dephases (e = I), A = £, corresponds 
to the classical (incoherent) case: G = X/z = (,/z. For 
strong disorder and weak coupling, A > Q, the trans- 
mission amplitude 7j becomes suppressed by localization 
effects. Then, T; consists of the product of the dephasing 
probability e and a factor £/a that counts the number 
of channels available for the interchain transport on the 
scale Q where the wavefunctions have appreciable am- 
plitude, multiplied by exp(— A/£) the ID coherent trans- 
mission amplitude between subsequent dephasing events, 
yielding G = e({/a)(X/z) exp(-A/C) = (C/z) exp(-A/C). 
Thus, for weak coupling, the classical expectation value 
for the conductance is suppressed by ID localization ef- 
fects. Delocalized carrier transport remains piecewise ID 
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FIG. 3: Metal-insulator phase-diagram of coupled ID chains 
as a function of interchain coupling and intrachain localiza- 
tion length. Below the drawn line the system is insulating, 
above it is metallic. Dots: disorder driven MIT in conducting 
polymers. The inset gives a graphical representation of the 
leitmotif behind the MIF criterion discussed in the text. 



and the length scale of interest is A with a ID transport 
time T\ = (C/vf) exp(A/£). At t\ scale, only fraction e of 
the carriers is involved in interchain transport and thus 
crossover to full 3D conduction occurs at T3D = T\/e. 

The above results demonstrate that in the presence of 
dephasing, the dc conductivity remains finite for z — > oo. 
This holds for any dephasing mechanism; in particular 
inelastic events due to finite temperature^ As argued 
above, "elastic dephasing" is a necessary condition for the 
formation of a truly metallic state. In our ID model, the 
interchain transfer couples a \1/id,i to a different 'J'id^ 
on lengthscale A and this gives a typical level splitting eA 
of these states. When the finite width of the ID energy 
levels (Thouless energy) on length scale A, given by h/r\, 
exceeds this level-splitting (bandwidth) energy-matched 
(elastic) interchain transitions are possible and hence a 
metallic phase can be formed, as indicated in the inset 
of Fig. [3J Applying the relation hvp ~ g^ 1 valid for 
ID conductors, the criterion for derealization translates 
into a/C = (B/E F ) 2 < eln(l/e). Alternatively, one can 
use the condition G > 1 for a metal on a characteristic 
length scale. When applied to the carrier fraction that is 
effectively involved in the interchain transport on length 
scale A this gives G e (X) — G(X)/e = cxp(— X/Q/e > 
1, leading to the same result. The phase diagram for 
coupled disordered chains as function of disorder (C/a) 
and interchain coupling (e) is shown in Fig. [3j 

Let us now apply our model to the metallic polymers. 
First we derive e using a tight-binding approximation. 
Interchain charge transfer would, in the absence of dis- 
order, lead to an interchain bandwidth oc However, 
as a result of disorder electronic states ^id on adjacent 
chains typically have an energy mismatch w A/2. Pro- 
vided that I± < A/2, the bandwidth available for inter- 
chain transport reduces to 4/^ /A which corresponds to 
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an interchain coupling e = 4(/_l/A) 2 governed by both 
interchain transfer and disorder. Next, we derive the rel- 
evant transport time (energy) scales. Crossover to full 
3D conduction occurs on time scale T3D = T\/e. Fol- 
lowing the reasoning for the MIT criterion, the band- 
width for interchain transport is eAexp(— A/£) and, us- 
ing tight-binding, the effective mass for 3D transport be- 
comes to* = 8?i 2 exp(A/C)/(eA6 2 ). ID localization sup- 
presses the transmission probability on scale A and the 
effective carrier density involved in 3D transport reduces 
to n e g = nexp(— A/£). This gives: 



1"3D 



(7rft/2JisC)exp ( a /4/ 2 ff 2 C 3 ), 



(2) 



2 2 6 2 / 2 g(/2h 2 £() exp (-a/4/ 2 5 2 C 3 ) . (3) 



For WT3D < 1 , the 3D transport channel is open and the 
conductivity is a(ui) = U3d/(1 + itoTsu)- At intermediate 
energy, ^/t 3 d < fouj < A, the ID signature of delocalized 
carrier transport becomes visible. Screening in ID con- 
ductors is poor and the dynamic response is described 
as cr(u>) — <r n j(cj)/(l + a n i(u>) / (iuiC)) , where a n i(uj) = 
oa/(1 + ioJT\) the non-interacting case and capacitance 
C s» e 2 g quantifies the ID screening efficiency^ For 
fiw > A the isolated chains' response (Lorentz oscilla- 
tor) is obtained: <j(lu) = oo/[l + i((u 2 — w 2 )/w) T o] with 
luq a measure for level separation. In case the time scales 
are sufficiently separated, the total optical response is 
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1 + icjT 3D 1 + iu)T\ 



iujC l + l — — -^Tq 



(4) 



Here 03D = £o^p T 3B, <J\ — <jzd because the factors e in t\ 
and bandwidth Aexp(— A/£) cancel and <7o = £o^p r o — 
ne 2 To/m the intrinsic intrachain conductivity. 

Here 03D = £ow 2 T3D, <ta = a 30 because the factors e 
in ta = eT3D and oj 2 A = Wp/e cancel and Co = £o^p T o = 
ne 2 To/m the intrinsic intrachain conductivity. 



Indeed, Fig. ^ demonstrates our theory gives an excel- 
lent description of the experimental results, using tq — 
3 fs, t\ = 0.4 ps, t 3D = 5 ps, 03D = o\ = 220 S/cm, 
cr = 380 S/cm, uj = 0.23 eV, and C = lO" 10 F/m. 
Thus, A/C = 3.5 and e = 0.08, which assuming a = 4 A 
leads to C = 1.4 nm and A = 5.0 nm. 

We conclude by addressing the scaling relation be- 
tween lo v and T3D- Eqs. and J3J express the ex- 
treme sensitivity of the 3D carrier dynamics for disorder 
and interchain coupling. For increasing disorder (lower 
C) and decreasing coupling (lower 7j_) the transition time 
to 3D conduction increases exponentially, while simulta- 
neously the plasma frequency decreases at the same rate. 
Let us consider the case of a disorder driven MIT. We 
chose To = 1.2-4.0 fs as typical range and ( — 4i>fTo, 10 
to reproduce the experimental results in Fig. [21 using 
known microscopic parameters for conducting polymers^ 
Ij_ = 50 meV, a = 4 A, b = 3 A, g = 0.5 eV -1 a~\ and 
n = 7 x 10 27 m" 3 , which corresponds to optimum doping 
of 0.25 carrier per monomer. A similar calculation can be 
performed for a MIT purely driven by interchain trans- 
fer, but agreement with experiment can be achieved only 
for an unrealistically large change in Jj_ of 20 to 90 meV. 
This agrees with the general consensus that the MIT in 
polymers is disorder driven. Using the above parameters, 
we have calculated the phase-diagram for the disorder- 
driven MIT, see Fig. [31 illustrating the transition to the 
metallic state. Indeed it shows the conducting polymers 
are on the edge of being metallic or insulating. 

In summary, we have discussed a transport model for 
weakly coupled disordered chains. Under appropriate 
conditions, (B / E-p) 2 < eln(l/e), a crossover from ID 
to 3D transport occurs and a metallic phase is formed. 
Expressions for the optical conductivity, 3D scattering 
time and plasma frequency have been derived. Based on 
established microscopic parameters, the model quantita- 
tively explains broad-band data on metallic PPy 
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